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Abstract 

Molecular dynamics simulations are carried out to study the translational and rotational diffusion 
of a single Janus particle immersed in a dense Lennard-Jones fluid. We consider a spherical particle 
with two hemispheres of different wettability. The analysis of the particle dynamics is based on the 
time-dependent orientation tensor, particle displacement, as well as the translational and angular 
velocity autocorrelation functions. It was found that both translational and rotational diffusion 
coefficients increase with decreasing surface energy at the nonwetting hemisphere, provided that 
the wettability of the other hemisphere remains unchanged. We also observed that in contrast 
to homogeneous particles, the nonwetting hemisphere of the Janus particle tends to rotate in the 
direction of the displacement vector during the rotational relaxation time. 

PACS numbers: 68.08.-p, 66.20.-d, 83.10.Rs 
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I. INTRODUCTION 


Understanding the processes of self- and directed-assembly of nanoparticles is important 
for design of various nanostructured materials with advantageous mechanical and optical 
properties jl] . Recent progress in synthesis of Janus and asymmetric patchy particles, that 


can form a variety o_ 
an industrial scale [h 


-5], 


makes their production feasible at 


predictable superstructures 

7]. The dynamics of assembly depends on the particle structure, dis¬ 
tribution of wettability and charges, and boundary conditions at the particle surface in a 
solvent. The hydrodynamic boundary conditions are usually specified via the so-called slip 
length, which is defined as an extrapolated distance with respect to the interface where the 
velocity profile vanishes js]. In general, it was shown that the slip length is a tensor quan¬ 
tity for flows over anisotropic textured surfaces [ 9 -12]. Transport properties of chemically 
homogeneous nanoparticles were recently studied in cylindrical pores filled with a wetting 


fluid 


13] and in an atomistic solvent confined between flat walls 14]. However, the diffusion 


process of individual Janus particles even in the absence of flow or confinement remains not 
fully understood. 

At the continuum level, the hydrodynamic behavior of Janus particles was studied in 


the low Reynolds number regime 


15 


than the particle size, the translational velocity o 
the moments of the force density on its surface 


18]. It was found that when the slip length is smaller 
' the spherical Janus particle is coupled to 
151 ] . Moreover, it was shown that due to 


its asymmetrical boundary conditions, Janus particles can migrate parallel to the velocity 


gradient in shear flows 


15]. In the limit of small slip lengths, the torque and force on 


spherical Janus particles with either discontinuous binary or continuously patterned surfaces 
were computed [16|,ll7]]. More recently, torques and forces on pill-shaped Janus particles with 
stick and slip boundary conditions were analyzed using the boundary integral method fis| . 
In particular, it was demonstrated that depending on the aspect ratio and slip length, the 
force transverse to the direction of a uniform flow can change sign 18| . In the present work, 
the diffusive motion of a Janus particle in a quiescent fluid was studied using molecular 
dynamics (MD) simulations. 

In recent years, Brownian motion of spherical particles with uniform wettability was 


extensively studied by MD simulations 


19-[26]. It was shown that a number of parameters, 


e.g., the mass and size ratios of a particle and a solvent molecule, interaction potentials, 
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boundary conditions and the system size, can influence the diffusion process. Notably, 



Furthermore, in the case of a smooth spherical particle, the Stokes-Einstein relation for the 

n 

diffusion coefficient holds for the system with slip boundary conditions [19[. However, if the 
particle consists of a cluster of atoms and their interaction energy with fluid molecules is 


sufficiently large, then stick boundary conditions apply at the particle surface 


M 


25|. 


In this paper, the results of molecular dynamics simulations are reported for the trans¬ 
lational and rotational diffusion of individual Janus nanoparticles in a quiescent solvent. In 
our study, Janus particles are modeled as spheres with two sides of different wettability, and 
their behavior is contrasted with that of uniformly wetting and nonwetting particles. We 
find that both translational and rotational diffusion coefficients are larger for Janus parti¬ 
cles with lower surface energy at the nonwetting hemisphere, provided that the wettability 
of the other hemisphere remains the same. It also will be shown that long-time algebraic 
decay of the velocity autocorrelation functions is recovered for uniform and anisotropic par¬ 
ticles. Interestingly, the diffusive motion of Janus particles involves a subtle correlation 
between translational and angular velocities; i.e., a non-zero velocity component parallel to 
the particle symmetry plane causes a finite rotation to reduce drag. 


The rest of the paper is structured as follows. The description of the molecular dynamics 
simulation model and parameter values is given in the next section. The analysis of the 
fluid density profiles, boundary conditions, particle displacement and rotation, as well as 
translational and angular velocity autocorrelation functions is presented in Sec. lIIIl The 
conclusions are provided in the last section. 


II. MOLECULAR DYNAMICS SIMULATION MODEL 


The diffusive motion of a Janus 
molecular dynamics simulations [32] 


particle in an explicit solvent was studied using classical 
In the explicit solvent model, the pairwise interaction 
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between fluid monomers is specified via the truncated Lennard-Jones (LJ) potential 


V LJ {r)= 4 e 


a 


L\r 


12 


- < 7) 


(1) 


where e is the depth of the potential well and cr is the separation distance at which the 
LJ potential is zero. The cutoff radius was fixed to r c = 2.5 a for all types of interactions. 
The system consists of 46 531 fluid monomers and a Janus particle, which are confined into 
a cubic box of linear size 39.62 a. All simulations were performed at a constant volume 
with the fluid density p = 0.75 cr -3 and temperature l-le/ks, where ks is the Boltzmann 
constant. The fluid temperature was maintained by the Nose-Hoover thermostat with the 
coupling parameter of 1.0. Periodic boundary conditions were applied in all three directions. 
The equations of motion were integrated numerically using the Verlet algorithm 30] with a 
time step A tuD = 0.005 r, where r = a^frnfe is the time scale of the LJ potential. 

The Janus particle was modeled as a convex polyhedron, the so-called Pentakis snub 


dodecahedron, which consists of 140 triangular faces, 210 edges, and 72 vertices 


31 ]. Out of 


the total 210 edges, there are 60 short (0.86 cr) and 150 long (0.99 a) edges. The radius of a 
sphere that contains all vertices is 2.08 cr. I 11 our model, 72 LJ atoms were fixed relative to 
each other at the vertices of the polyhedron; and, during the simulations, the Janus particle 
was moving as a rigid body under forces from the surrounding fluid. The particle atoms 
do not interact with each other and their size is the same as the size of fluid monomers. 
The particle mass is M = 50.0 m. The surface of the Janus particle was divided into two 
hemispheres of different wettability, as shown in Fig.[0 The LJ interaction energy between 
fluid monomers and particle atoms at the wetting hemisphere was fixed to e p f = 1.0 e, and 
at the nonwetting hemisphere £ p f was varied in the range 0.1 e ^ £ p f ^ 1.0 e. The two 
limiting cases of uniformly wetting and nonwetting particles correspond to the values of the 
interaction energy with the fluid molecules of 1.0 £ and 0.1 £, respectively. 


The Janus particle and the surrounding fluid phase were first equilibrated for about 
2 x 10' MD steps (or 10 5 r) and then the data were accumulated during the next 6 x 10 8 
MD steps (or 3 x 10 6 r). During the production runs, the position and orientation of the 
Janus particle, the velocity of all particle atoms as well as the position of the center of mass 
of the solvent were saved every 20 MD steps. To accumulate statistics, the simulations were 
repeated in forty independent systems for each value of e p f. In each case, the velocities of 
fluid monomers were initialized using the Maxwell-Boltzmann probability distribution with 
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zero total linear and angular momenta 32], We also checked that after equilibration, the 
Maxwell-Boltzmann distribution for the translational and angular velocities of the Janus 
particle is recovered. 


III. RESULTS 


It has long been recognized that in the presence of a smoot 
tend to form several distinct layers parallel to the substrate 


i substrate, the fluid monomers 


33]. The amplitude of the fluid 


density profiles gradually decays away from the substrate to a uniform bulk density profile. 
Similarly, in the case of a spherical particle, the spatial distribution of fluid monomers 
around the particle consists of several concentric layers 38]. In our simulations, the radial 
distribution function for fluid monomers is shown in Fig.[2]for uniformly wetting (e p f = 1.0 e) 
and uniformly nonwetting (e p f = 0.1 e) particles. In both cases, the fluid density profiles 
decay to a uniform profile within about (3 — 4 ) a from the particle surface. As expected, 
fluid layering is more pronounced for the uniformly wetting particle. In simulations with a 
Janus particle, a reduced value of the surface energy at the nonwetting hemisphere resulted 
in a lower amplitude of the local fluid density profile above the nonwetting side, while the 
density oscillations remain large above the wetting hemisphere where the interaction energy 
is fixed to e (see the inset in Fig. [2]). 

The flow boundary conditions at the particle surface depend on the atom density and 
the LJ interaction energy between particle atoms and fluid monomers. In order to estimate 
the Navier slip length, we performed a set of separate simulations of fluid flow (without a 
particle) confined between flat crystalline walls with the same density as the Janus particle. 
In this setup, the fluid phase with density p = 0.75 cr ~ 3 is in contact with (111) planes of 
the face-centered cubic (fee) lattice. Similar to the previous studies 39s,|40j, the Poiseuillc 


flow was induced by a constant force f x = 0.005 e/o applied to each fluid monomer in the 
direction parallel to the walls. The slip lengths were computed from a parabolic fit of the 
averaged velocity profiles for different values of the wall-fluid interaction energy. The results 
are summarized in Tabled] It is evident that at high interaction energies, the slip length 
is about the particle radius, while at low energies the slip length becomes larger than the 
particle diameter. Note also that for flows over curved surfaces the slip length should be 


modified by the local radius of curvature 


34 


31 ] 
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In a quiescent solvent, the Janus particle undergoes translational and rotational diffusion 
under random forces from the fluid monomers. We start the discussion of our results by 
considering the position vector of the particle center of mass. The mean square displace¬ 
ment (MSD) for uniformly (non)wetting and Janus particles is plotted in Fig. E] for different 
surface energies. In all cases, the center of mass of the fluid phase was subtracted from the 
particle position. It can be clearly seen that after the ballistic regime, the MSD increases 
linearly with time (as indicated by the dashed line in Fig. [3]). As expected, the slowest and 
fastest diffusion correspond to uniformly wetting and nonwetting particles, respectively. In 
between, the displacement increases with decreasing surface energy at the nonwetting hemi¬ 
sphere of Janus particles. Overall, this result is not surprising and can be expected for any 
spatial distribution of binary wettability on the particle surface. The translational diffusion 
coefficient was computed from the relation r 2 = 6 D t t at times t > 10 3 r (shown in the upper 
inset of Fig.[3j). 

The rotational dynamics of nanoparticles can be probed by analyzing the decay of the 
autocorrelation function of the unit vectors e* that are fixed relative to the particle. In 
what follows, we choose the vector ei along the main axis of symmetry of the particle (see 
Fig. CD)- The autocorrelation function (e^O) • e^f)) is plotted in Fig. [4] for the Janus and 
uniformly wetting and nonwetting particles. One can observe a similar trend as with the 
translational diffusion; namely, with increasing surface energy at the nonwetting hemisphere, 
the rotational relaxation slows down, and the two limiting cases of uniformly (non)wetting 
particles correspond to (fastest) slowest decay. A typical rotational relaxation time, r r , was 
estimated from an fit to a simple exponential decay, i.e., 

(ei(0) • ei(t)) = exp(—t/r r ). (2) 

The results are reported in the inset of Fig.[4l It can be observed that r r gradually increases 
with increasing surface energy e p f, which implies that the rotational diffusion coefficient is 
larger for Janus particles with lower surface energy at the nonwetting hemisphere. 

In order to emphasize the difference between uniform and anisotropic particles, we next 
consider the translational and angular velocity autocorrelation functions (VACF) defined 
as (v(0) • v(t)) and (u>(0) • cj(t)), respectively. First, we plot the translational VACF in 
Fig. [5] for uniformly (non)wetting and Janus particles. Note that in all cases, the VACF is 
determined by the thermal velocity v 2 = 3 fegT/M at short times, while the algebraic decay 
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results 


e si 

ope 

-1 

22, 

25, 

27 


1.5 is recovered at long times, which is consistent with the previous MD 


27]. In the intermediate regime, the VACF is smaller for Janus particles 


with higher interaction energy at the nonwetting hemisphere. Similarly, the angular VACF 
for Janus particles gradually varies between the two limiting cases of uniformly wetting and 
nonwetting particles (shown in Fig.E]). Although the statistics at long times is limited, it 
appears that the angular VACFs exhibit a power-law decay with the slope —2.5, in agreement 


with earlier predictions 


28 


29]. 


The autocorrelation functions discussed above do not necessarily reflect the broken sym¬ 
metry of Janus particles; in principle, one can expect an enhanced diffusion with decreasing 
surface wettability regardless of its spatial distribution. In what follows, we explore a corre¬ 
lation between translational and rotational motion of Janus particles. It is intuitively clear 
that if a constant force is applied to a Janus particle in a viscous fluid then the nonwet¬ 
ting face will be oriented along the direction of motion in order to reduce drag. In analogy, 
when a Janus particle moves under random forces from the surrounding fluid, its nonwetting 
hemisphere tends to rotate in the direction of the displacement vector. This effect can be 
quantified by computing a scalar product of the displacement vector of the center of mass 
during the time interval t and a vector difference between unit vectors ei(t) and e^O), and 
then plotting it as a function of time (not shown). However, we found that it is difficult to 
obtain good statistics for this quantity and, more importantly, to provide a clear physical 
interpretation of the effect. 


Alternatively, a more direct way to probe the correlation between translation and orien¬ 
tation of Janus particles is to compute a rotation angle along the direction of displacement. 
More specifically, we first choose a time interval t, compute the displacement vector of the 
center of mass Ar, determine angles between (ei(0), Ar) and (ei(t), Ar), and then plot 
the difference between these angles as a function of time. In other words, we measure how 
the unit vector ei rotates along but not around the displacement vector. The averaged 
rotation angle is plotted in Fig.|T]as a function of time for the Janus and uniformly wetting 
and nonwetting particles. It can be clearly observed that, within statistical uncertainty, 
particles with uniform distribution of wettability do not, on average, rotate along the dis¬ 
placement vector. In contrast, Janus particles undergo a finite rotation, which increases 
with decreasing surface energy at the nonwetting hemisphere. In each case, the rotation 
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angle is maximum at times of about the rotational relaxation time (see inset in Fig.[4j). Note 
also that the maximum rotation angle in Fig. [7] occurs at larger times for Janus particles 
with higher surface energy, which correlates with a slower relaxation of the autocorrelation 
function (ei(0) • ei(f)) for Janus particles with higher e pf reported in Fig.[H As expected, 
at longer times the correlation gradually decays, although it remains nonzero even at times 
of about 10 4 r. 

To further quantify the effect of rotation along the displacement, we also measured a 
correlation between the translational and angular velocities. The corresponding correlation 
functions were defined as a product (vi(O)-ojj(t)), where the velocity components with indices 
i,j = 1,2,3 were computed along the unit vectors e* (see Fig.[T]). As a representative case, 
we consider a Janus particle with the largest difference in wettability, (1.0e, 0.1 £), when 
the correlation effect is most pronounced. The averaged correlation functions are presented 
in Fig. [8] for selected pairs of indices. It can be seen that the only non-zero correlation 
functions are (^(O) ■ u}s(t)) and (^(O) Physically, it means that if a Janus particle 

has a translational velocity component in a direction parallel to the symmetry plane then 
it will tend to rotate its nonwetting hemisphere towards the displacement. Finally, we note 
that for uniformly wetting/nonwetting particles there is no preferred rotation with respect to 
translation, and, therefore, the cross correlation functions are zero at all times (not shown). 


IV. CONCLUSIONS 

In summary, we investigated the translational and rotational diffusion of individual Janus 
nanoparticles in a dense fluid using molecular dynamics simulations. In our model, the 
Janus particle consists of atoms arranged on a sphere, which is divided in two hemispheres 
of different wettability. The interaction energy between particle atoms and fluid molecules 
is fixed at the wetting hemisphere, while the wettability of the other side is a variable 
parameter. The diffusion process of such particles is compared against two limiting cases of 
uniformly wetting and nonwetting particles. 

We found that the diffusion coefficient of Janus particles is bounded between the two 
limiting cases and it increases with decreasing surface energy at the nonwetting hemisphere. 
The analysis of the translational and angular velocity autocorrelation functions showed that 
the exponents of the long-time power-law decay are the same for Janus and uniform particles. 
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Perhaps most interestingly, we observed a correlation between translation and rotation of 
Janus particles; in other words, when a particle spontaneously moves in a particular direction 
it tends to reduce drag by rotating its nonwetting side towards the displacement vector. This 
effect is absent for homogeneous particles. 

In the future, it will be instructive to extend these results to a system of many particles 
and study in detail the dynamics of self-assembly of Janus particles into complex structures 
in the bulk or at liquid-liquid interfaces. 
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Figures 
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FIG. 1: (Color online) A snapshot of the Janus particle with wetting (blue atoms) and nonwetting 
(green atoms) hemispheres. The particle consists of 72 atoms rigidly fixed on vertices of a convex 
polyhedron (see text for details). Note that particle atoms are not drawn to scale. The unit vectors 
ei, e 2 and e 3 form an orthonormal basis in a reference frame of the particle. 
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r/a 

FIG. 2: (Color online) Averaged fluid density profiles as a function of the distance from the particle 
center for uniformly wetting (o) and uniformly nonwetting (v) particles. The LJ interaction energy 
between particle atoms and fluid monomers is £ p f = 1.0£ (o) and e p f = 0.1 e (V). The inset shows 
fluid density profiles near the wetting hemisphere £ p f = 1.0 e (x) and nonwetting hemisphere 
e p f = 0.1 e (o) of the Janus particle. 


13 


















FIG. 3: (Color online) The mean square displacement of particles as a function of time for the 
surface energies at both hemispheres from top to bottom (0.1 s, 0.1 s), (1.0s, 0.1s), (1.0 s, 0.3 s), 
(1.0 s, 0.5s), (1.0s, 0.7s), and (1.0s, 1.0s). The dashed line with unit slope is plotted for reference. 
The lower inset shows an enlarged view of the same data at large times. The upper inset shows 
the diffusion coefficient as a function of s p f, while the other hemisphere is wetting, i.e., s p f = 1.0 s. 
The case s p f = 0 corresponds to a uniformly nonwetting particle (0.1s, 0.1s). 
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FIG. 4: (Color online) The autocorrelation function (ei(0) • ej (t)) for the surface energies at both 
hemispheres (O.le, O.le), (l.Oe, O.le), (l.Oe, 0.3e), (l.Oe, 0.5e), (l.Oe, 0.7e), and (l.Oe, l.Oe) 
from left to right. The inset shows the rotational relaxation time r r (in units of r) as a function of 
the surface energy at the nonwetting hemisphere. 
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FIG. 5: (Color online) The velocity autocorrelation function (v(O)-v(t)) for the interaction energies 
at both hemispheres (0.1 e, 0.1 e), (l.Oe, 0.1 e), (l.Oe, 0.3e), (l.Oe, 0.5e), (l.Oe, 0.7e), and (l.Oe, 
1.0 e) from top to bottom. The straight line with a slope —1.5 is shown as a reference. 
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FIG. 6: (Color online) The angular velocity autocorrelation function (u>(0) • u>(i)) for the LJ 
interaction energies at both hemispheres (0.1 e, 0.1 e), (l.Oe, 0.1 e), (1.0s, 0.3e), (l.Oe, 0.5e), 
(l.Oe, 0.7e), and (l.Oe, l.Oe) from top to bottom. The black line with a slope —2.5 is plotted for 
reference. 
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FIG. 7: (Color online) The rotation angle along the displacement vector of the center of mass for 
the surface energies at both hemispheres (1.0e, 0.1 e), (1.0e, 0.3e), (1.0 e, 0.5e), (1.0s, 0.7e) from 
top to bottom. The data for uniformly wetting (1.0 e, 1.0 e) and nonwetting (0.1 e, 0.1 e) particles 
are indicated by the black and red curves, respectively. 
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FIG. 8: (Color online) The translational times angular velocity correlation functions (u*(0) • u>j(t )) 
for the indicated pairs of indices. The velocity components are computed along the unit vectors e* 
(see Fig.|T]). The data are reported for a Janus particle with the LJ interaction energies at both 
hemispheres (1.0e, 0.1 e). 
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TABLE I: The slip length L s /a computed in a Poiseuille flow confined by flat crystalline walls with 
density 1.63 for the tabulated wall-fluid interaction energies. The error bars for the slip length 
are about 0.5 a. 


£pf/ £ 

0.1 

0.3 

0.5 

0.7 

1.0 

L s /a 

9.0 

6.4 

4.6 

3.4 

2.0 
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